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METHOD FOR CONTROIjLING COMBUSTION IN AN INTERNAL 
COMBUSTION ENGINE AND PREDICTING PERFORMANCE AND 
EMISSIONS 

Abstract 

This disclosure teaches a method of controlling a 
direct injection internal combustion engine and 
predicting the behavioxir of a direct injection 
internal combustion engine. An estimation of 
initial cylinder press\are, air flow and EGR flow 
(if applicable) is used to establish a system 
that provides engine behaviour by integrating an 
injection module, combustion module and engine 
control module to provide data indicative of 
engine behaviour such as brake torque and power, 
air flow, EGR flow, cylinder pressure, brake 
specific fuel consunqption, start of combustion, 
heat release rate, ttirbo- charger speed and other 
variables. These values can then be used to 
adjust commanded variables such as start of 
injection, commanded pulse width, rail pressure 
to meet operator demand. Also the output data 
can be used as a tool to determine how a 
conceptualised engine design will behave. This is 
particularly useful for gaseous- fuelled internal 
combustion engines where cylinder pressure 
influences behaviour of injected gases in light 
of the fact that rail pressure and cylinder 
pressure are, generally, of a similar magnitude. 
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METHOD FOR CONTROLLING COMBUSTION IN AN INTBRNiUj 
COMBUSTION ENGINE MID PREDICTING PERFORM2VNCE AND 

EMISSIONS 

Field of the Inventlcn 

The present Invention relates to a method 
for controlling coitibustion, performance and 
emissions of a direct injection internal 
5 combustion engine and predicting the behavior of 
a direct injection internal combustion engine. 
An embodiment of the subject method is 
particularly useful for gaseoxxs- fuel led direct 
injection internal combustion engines. 

10 

Background of the Invention 

Accurate prediction of the behavior of 
internal combustion engines without engine 
testing can be valuctble. First, predicting how 

15 an engine will perform under a set of conditions 
can reveal valuable information without having to 
incur the cost and time involved in testing an 
actual engine. Second, e^qpensive prototype 
refinement is not unusual when testing a new 

20 engine design. Predicting behavior of a concept 
without having to manage prototype development 
can be a helpful step in engine research 
accommodated by methods of predicting engine 
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helps overcome these problems by initially 
testing new ideas for a new engine design or 
simulating an existing engine architecture prior 
to any investment in costly prototype parts or 
5 time consuming test cell work. Finally, methods 
of predicting engine behavior can allow for 
accelerated testing of a potential system without 
limitation of manpower and facility availability. 
While not a substitute for actual testing, 

io accurate prediction can eliminate many otherwise ' 
reasonable ideas before an expensive test engine 
is built or before valuable test cell time is 
used to test an existing engine. 

Additionally, as computer processing power 

15 continues to inprove, open loop controls using 
engine prediction data can eliminate both closed 
loop control methods and reliance on expensive 
monitoring equipment and the associated costly 
electrical systems, both of which are prone to 

20 wear. 

Finally, as emissions controls become 
increasingly important for environmental reasons, 
it is valuable to provide a controller for closed 
loop on-board control adjusting emissions in 

25 response to predicted behavior for a set of user 
demanded operational conditions. Also, as is the 
case in engine design noted above, a predictive 
system can help to filter out new designs that 
fail to meet acceptable emission levels. 

3Q Prior art in this area has provided, in 

general, systems that predict aspects of engine 
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operation but fail to satisfactorily integrate 
these to provide a predictive engine. By way of 
example, combustion models have been developed 
that predict combustion within a combustion 
chamber on the basis of operator inputted values. 
The same is true of emissions models, injector 
models and engine cycle models. 

Examples of engine combustion models 
development are summarized in Heywood, 
Jntemal Coznbustion Engine Fundamentals (New 
York: McGraw-Hill, 1988) • The combustion model 
is typically integrated into the engine cycle 
model. Often accompanying these combustion models 
are emission formation models that provide an 
indication of NOx and particulate matter 
emission. Engine cycle simulation is used to 
compute thermodynamics and gas dynamics in the 
engine. GT- Power™ is an example of an engine 
cycle process modeling module. For engine cycle 
prediction, several methods have been provided in 
such systems as StarCD™, Fluent™, FIRE™ and 
VECTIS™. Hydsim^"^ software is an example of a 
prior art injector simulation model* The software 
can compute injection delay, injection rate shape 
by modeling details of the injector, including 
its hydraulic system, mechanical system and 
electronic system. 

Each of these components has been developed, 
to a greater or lesser extent, but, in general, 
remain independent models of one aspect of engine 
behavior. By way of example, consider A. Chow et 
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al.. Thermodynamic modeling of complete engine 
systems - a review (London: Professional 
Engineering Piiblishing, 1999) and T. Morel et 
al./ ^Virtual Bngine/Powertreiin/Vehicle'' 
5 Siinulation Tool Solves Complex Interacting System 
Issues (Warrendale, PA: Society of Automotive 
Engineers International, 2003) 2003-01-0372. In 
each case, examples of models used to describe 
aspects of an engine are provided. 
10 Chow discusses the system integration of 

models which determine engine behavior for a 
specified rate of combustion. These raodels 
quantitatively represent heat transfer, fluid 
flow, t\irbo-machinery, and emissions production. 
15 AS one example of a system integration, Morel 

describes software which, for a specified rate of 
combustion, provides a system integration of 
separate models for heat transfer, fluid flow, 
turbo -machinery, and hydraulic, electrical, 
20 thermal, mechanical r and controls elements of the 
engine system. Other prediction methods 
discussed by Morel include components with 
turbochcuiHer, supercharger and intercooler 
support; co-simulation options; transient pipe 
25 teitqperature analysis; variable pipe length; 

mult i -user database and dynamic valve motion. The 
methods taught emphasize the relationship of the 
engine valves with other components. 

Generally, the integrations discussed do not 
30 explicitly include the interaction between fuel 
injection behavior and the rate of combustion. 
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Full engine simulation integration benefits if 
the method allows for the interdependence of 
injector performance, combustion pressure and 
combustion rate. 
5 Generally, limited development in integrated 

engine predictive methods for direct injection 
engines has been due, in part, to individual 
component modeling for direct injection diesel 
fuelled engines that has been satisfactory for 

10 the piorposes of modeling such engines. 

For example, the changing conditions inside the 
combustion chamber caused by combustion in a 
diesel engine has limited influence on the 
injection of diesel fuel due, largely, to the 

15 extremely high pressures at which diesel fuel is 
injected, so an integrated model that adjusts the 
fuel deliveary to a cylinder in response to the 
cylinder pressure does not suffer to any 
significant extent by ignoring combustion of the 

20 fuel in the cylinder. That is, compared to the 
pressures established in the combustion chamber 
as a result of combustion events, the pressure at 
which diesel fuel is injected is many times 
higher and, therefore, less influenced by changes 

25 in combustion chamber pressure. 

However, ignoring this part of an engine 
system for a predictive model can be a problem 
when a gaseous fuel is used to operate high 
cott^ression ratio direct injection engines. For 

30 diesel engines operating on gaseous fuels, the 
fuel injection rate depends in general on 
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pressxzre in the combustion chaniber; however the 
rate of change of cylinder pressure depends on 
the injection rate, as well as on the turbo- 
machineiry and heat transfer components. Thus for 
5 a full system integration separate models for 
predicting injector behavior and combustion 
development are important. These models exist 
independently but have not been adequately 
integrated into a full system applicable to 

10 diesel engines, or engines with diesel engine 
compression ratios, where those engines operate 
on gaseous fuels. 

Most component models simulate aspects of an 
engine removed from the operation of the engine 

15 system as a whole. The interface between the 

modeled components and the rest of the system is 
based on engine test measurements or data. 
Therefore, the model is not a self -supported 
full -predictive model. The models are applied to 

20 understand the engine process for diagnostic 
purpose • 

As noted above, in injector computation, the 
model upstream conditions, such as rail pressure 
and temperature, are from measurements that, in 

25 general, do not consider fluctuations in 

pressure, which might be quite significant in 
some cases. This is especially important when gas 
fuel injection is involved. In combustion 
models, the intake conditions are usually from 

30 test measurements on a functioning engine- There 
is little consideration of the interaction 
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between the injection and the combustion. 
Likewise, in engine cycle simulation, the 
combustion process is often derived from the heat 
release rate measured on an actual engine. This 
5 ignores the effect of the intake conditions on 
the combustion process and simplifies the 
interaction between injection, combustion and air 
intedce axid management processes. It is believed 
that the interaction between components in the 

10 engine system is very important and needs to be 
well understood in order to shorten engine 
development phase. However, it is also very 
complex. Most engine models do not reflect the 
complex interactions in engine system. Models 

15 are detached from the system reducing a model's 
predictive capability. 

Moreover, until recently computer power, has 
not been adequate to accurately integrate 
modeling coirq?onents to predict overall engine 

20 behavior in a time frame useful for the purposes 
contetrplated. Developments in computer 
processing power have provided the necessary 
computational capacity to model an integrated 
engine • 

25 Efforts to integrate sub-models together 

have, in general, not resulted in integrated fuel 
systems, combustion system and engine breathing 
systems (as well as emissions systems) - In 
particular, modeling efforts have not been 

30 developed for such a predictive integrated system 
for a gaseous- fuelled engine. Also, while the 
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focus of efforts to create engine predictive 
methods has been directed to the purpose of 
providing an engine diagnostics tool and 
predicting engine behavior, little effort has 
been made on the use of such methods for engine 
controls . 

The present method discloses an integrated 
engine prediction tool for controlling and 
predicting engine performance and emissions. 

Simmary of the Invention 

The present invention provides for an 
integrated engine model and method for predicting 
and controlling a direct injection fuelled engine 
and especially a gaseous -fuel led direct injection 
engine . 

A ftirther embodiment of the present 
invention includes A method is disclosed for 
controlling a direct injection internal 
combustion engine by calculating output data from 
system data. The method con^rises using system 
data which includes an injector geometry defined 
by a fuel Injector in fluid communication with a 
combustion chamber of the engine, an engine 
geometzy defined by the engine, an engine 
operating command based on a desired engine 
output, that is indicative of a commamded torque 
for the engine, an initial injector command, and 
an engine speed. The output data comprises an 
engine out value indicative of a brake torque 
delivered by the engine, and sua output control 
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parameter indicative of a cylinder pressure 
value* When the output data does not satisfy a 
predetermined relationship, it is recalculated 
with a determined subsequent injector command. 
This predetermined relationship compaxes the 
output data, the engine operating command and a 
demanded control parameter. The demanded control 
parameter is indicative of a maximum cylinder 
pressure value. VOien the output data satisfies 
the predetermined relationship, the injector is 
commanded to provide fuel to the combustion 
chamber according to a last injector command. 

In a further embodiment of method, the 
system data can also comprise an ambient 
temperature, an ambient pressure and an engine 
speed, which in a preferred embodiment are 
measured prior to calculating the output data. 
In a further embodiment of the present method, 
the engine operating command is converted to the 
commanded torque based on an operator selected 
pedal position. 

In a further embodiment of method, the 
demanded control parameter is determined from the 
engine operating command. As well, the output 
control parameter is indicative of calculated 
emissions and the demanded control parameter is 
indicative of a demanded emission. The calculated 
emissions and the demanded emissions are 
indicative of NOx and/or particulate matter 
concentration. 

A further embodiment of method discloses the 
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predetermined relationship comprising comparing 
the output data with a previously calculated fuel 
consumption determined from the system data using 
the initial or subsequent injector commands. The 
output data further comprises the calculated fuel 
consumption. To meet the predetermined 
relationship, the calculated fuel consumption is 
less than the previously calculated fuel 
consiimption. 

In a further embodiment of method, the 
initial injector command comprises commanding 
initial start of injection, an initial pulse 
width and an initial rail pressure. The 
subsecpient injector command coii^rises commanding 
a subsequent start of injection, a subsequent 
pulse width and subsequent rail pressure. And 
the last injector command comprises commanding a 
last start of injection, a last ptilse width and a 
last rail pressure. 

A further embodiment of method discloses a 
variable geometry turbine. The system data may 
conprise variable turbine geometry positions. 
The injector commands Include the last variable 
geometry turbine positions when the output data 
satisfies the predetermined relationship. When 
the output data does not satisfy the 
predetermined relationship ^ a subsequent variedale 
geometry turbine position, with the subsequent 
injector command, is used to recalculating with 
the output data. 

In a further embodiment of the method the 
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engine comprises an exhaust gas recirculation 
(EGR) system. The system data further comprises 
an initial effective EGR valve flow area^ further 
comprising/ when the output data does not satisfy 
5 the predetermined relationship, determining a 

subsequent effective EGR valve flow area with the 
subsequent effective EGR valve flow area and the 
subsequent injector command, The output data 
being recalculated when the output data satisfies 

10 the predetermined relationship resulting in 
commanding the EGR system according to a last 
effective EGR valve flow area. 

In a further embodiment of the method, the 
output data is calculated using an injector 

15 module for calculating IM/CM data and IM/ECM data 
from an estimated initial cylinder pressure and 
IM system data. The IM system data is a subset 
of the system data. Also, a combustion module 
for calculating CM/ECM data from an estimated 

20 initial intake flow rate, the IM/CM data and CM 
system data is used. The CM system data is a 
subset of the system data. Also, an engine cycle 
module for calculating the output data from the 
CM/ECM data, the IM/ECM data and the ECM system 

25 data is used. The ECM system data is a siabset of 
the system data. Prior to calculating the output 
data, a small loop between the combustion model 
and the engine cycle module is used to provide a 
converged intake flow rate. A big loop between 

30 the injection, combustion and engine cycle 

modules is used to provide a converged cylinder 
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pressure, which then provides the cylinder 
pressure value. 

In a further embodiment of the present 
method, the IM/CM data can comprise an IM rate of 
injection and an IM start of injection. An IM 
fuel flow and an IM start of injection can also 
be used. As well, the CM/ECM data can comprise a 
CM heat release rate and a CM start of 
combustion. 

The output data is further calculated using 
an emissions module for calculating emissions 
data from ECM/EM data, the IM/BCM data and the 
system data. The output control parameter is 
indicative of calculated emissions and the 
demanded control parameter is indicative of a 
demanded emission. 

Further embodiments of the present method 
include the IM system data comprising an ambient 
pressure, the injector geometry, the engine speed 
and the injector command. The CM system data 
can further comprise the engine speed and the 
engine geometry. The ECM system data comprises an 
ambient pressure, an ambient temperature, the 
engine geometry and the engine speed. 

The method can be practiced wherein the 
engine is at least partially fuelled by a gaseous 
fuel and in a preferred embodiment the gaseous 
fuel is natural gas. The gaseous fuel can also 
comprise methane, ethanol and propane. In a 
further preferred embodiment the fuel is 
hydrogen. 
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A Ttiethod is disclosed for predicting output 
data generated by a direct injection internal 
coTttoustion engine, comprising at least one of 
engine torque, fuel consumption, power and 
5 cylinder pressure. The method comprises 

selecting system data which comprises an injector 
geometry indicative of a selected fuel injector, 
an engine geometry indicative of the engine, an 
engine speed and an injector command. The method 
10 also comprises estimating initial data, 

comprising an initial cylinder pressure and an 
initial intake flow. The output data is 
calculated using the system data and a converged 
cylinder pressiire and a converged intake flow. 
15 The initial data is used to provide the converged 
cylinder preesxire from the initial cylinder 
pressure, and provide the converged intake flow 
from the initial intake flow. 

In a further embodiment, the system data can 
20 comprise either an ambient pressure and/or an 
ambient teraperat\ire while the output data can 
comprise NOx and/or particulate matter 
concentration. As well, the injector command may 
comprise a commanded rail pressure, a commanded 
25 start of injection aawi a commanded pulse width. 
In a further embodiment of the present 
method, system data can coiiqprise a commanded 
variable geometry turbine position and an 
effective EGR valve flow area. When an effective 
30 EGR valve flow area is included, the system data 
fxirther cott^jrises an initial EGR flow. 
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In a further embodiment the output data is 
calculated using an injector module for 
calculating IM/CM data and IM/ECM data from the 
initial cylinder pressure and IM system data. 
5 The IM system data is a siabset of the system 
data. Also, a combustion module is used for 
calculating CM/ECM data from the initial intake 
flow, the IM/CM data and CM system data. The CM 
system data is a subset of the system data. Also 
10 an engine cycle module is used for calculating 

the output data from the CM/ECM data, the IM/ECM 
data and ECM system data. The ECM system data is 
a subset of the system data. With this, a small 
loop between the combustion model and the engine 
15 cycle module is used to converge to a converged 
intake flow. As well, a big loop between the 
injection, combustion and the engine cycle 
modules is used to converge to the converged 
cylinder press\ire. 
20 In a further embodiment of the method, the 

output data is further calculated using an 
emissions module for calculating emissions data 
from ECM/EM data, the IM/ECM data and EM system 
data. The EM system data is a subset of the 
25 system data. The output data cottqprises the 

emissions data, which, in a preferred embodiment, 
is indicative of NOx concentration and/or a 
pcurticulate matter concentration. 

In a further eitibodiment, the ECM/EM data 
30 comprises the converged cylinder pressure, the 
converged intake flow, an intake valve closing 
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time cylinder pressure and an intake valve 
closing time cylinder temperature. The intake 
valve closing time cylinder pressure and 
temperature are provided by the engine cycle 
module . 

In a further embodiment of the subject 
invention, the IM/CM data comprises an IM rate of 
injection and an IM start of injection. The 
IM/ECM data comprises an IM fuel flow and an IM 
start of injection. The CM/ECM data comprises a 
CM heat release rate and a CM start of 
combustion. The IM system data can comprise the 
injector geometry, the engine speed, the ambient 
pressure and the injector command • The CM system 
data comprises the engine speed and the engine 
geometry- The ECM system data comprises the 
engine geometry, the ambient pressure, the 
ambient temperature and the engine speed. 

The method can be practiced wherein the 
engine is at least partially fuelled by a gaseous 
fuel and, in a preferred example, the gaseous 
fuel is natural gas. The gaseous fuel can also 
comprise methane, ethanol auid propsuae. In a 
further preferred example the gaseous fuel is 
hydrogen . 

Brief Description of the Drawings 

FIG. 1 shows a flow diagram of an integrated 
engine model for use in engine control or 
predicting engine behavior, 

FIG. 2 shows a flow diagram demonstrating a 
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small iteration loop. 

FIG. 3 shows a flow diagram demonstrating a 

big iteration loop- 

PIG. 4 shows the data flow path of the 
initial big loop as shown in FIG. 3. 

PIG. 5 shows the data flow path(s) for the 
subsequent big loops subsequent to the initial 

big loop of FIG. 4, 

FIG. 6 shows the data flow path wherein the 
integrated engine model is used for engine 
controls . 



15 



r.^*.^A^^ Daacripti.oa «^ Preferred Bmbodunentisl 

In the present disclosure, a method is 
taught to help with engine control, particularly 
for gaseous- fuelled direct injection internal 
combustion engines. The method predicts internal 
combustion engine behavior. 

Referring to Fig. l, the flow diagram 
20 illustrates integrated engine method (IBM) 8. 
The flow diagram shows the integration of 
injector module (IM) 10, combustion module (CM) 
12, engine cycle module (ECM) 14 and emissions 
module (BM) 16, each of 'which calculates output 
25 data based on provided input data where the 

output data is descriptive of the behavior of the 
injector, the selected cojnbustlon mode for the 
chosen fuel, the engine behavior and the 
emissions generated by the engine of interest, 
system data 18 includes geometry parameters and 
operation parameters, to be described in greater 



30 
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detail below. System data 18 is used to initiate 
operation of lEM 8. System data 18 lists all 
data used for IBM 8. Subsets of system data 18, 
IM system data 20b, CM system data 21b, ECM 
system data 31 and EM system data 39, each 
provide system data 18 used for each module. 

As well, initial data, namely, initial 
cylinder pressure (P^) data 19 and initial intake 
air flowrate (Q^) . Exhaust gas recirciilation 
flow rate (g^) is also provided as initial data 
in those cases where exhaust gas recirculation 
(EGR) is used or considered. and are, 

collectively, CM input data 21 where EGR is used 
and fife alone where it is not used. 

Lines with arrowheads represent data paths 
for routing data throughout lEM 8 between modules 
10, 12, 14 and 16, as well as routing system data 
18, P^i data 19 and CM input data 21 into IBM 8. 
Data paths 20, 22 and 24b direct IM system data 
20b, data 19, data 23, respectively, to IM 
10. Data paths 26 and 28 direct IM/ECM/EM data 25 
and IM/CM data 27, respectively, from IM 10 as 
indicated. IM/BCM/EM data 25 is directed to BCM 
14 and EM 16. IM/CM data 27 is directed to CM 
12. IM system data 20b, a subset of system data 
18 con5)rises, injector geometry data {GEO^) , 
injector commands, in the embodiment shown rail 
pressure iP,^) , commanded start of combustion 
{SOJg^) , and commanded pulse width iPW^^) , as 
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well as and engine speed {RPM„,) and ainbient 
pressure of the engine of interest </*,„^) • 

Data path 24b originates from path 24 and 
directs P;, data 23 to IM 10 after con^jletion of 
an iteration of big loop as will be discussed 
below. 

Along with IM/CM data 27 and CM input data 
21 (directed through path 40) , as noted above, CM 
12 also takes CM system data 21b throu;gh path 30c 
and generates CM/ECM data 29. After an iteration 
through ECM 14, as will be described below, CM 12 
also uses ECM/CM data 35 directed through path 34 
replacing CM input data 21 on subsequent 
iterations . 

ECM 14 uses IM/ECM data 25, CM/ECM data 29 
and ECM system data 31 from paths 2Sb, 32 and 
30a, respectfully, to provide P^, data 23, ECM/EM 
data 33 and ECM/CM data 35 directed to IM 10, EM 
16, and CM 12 through paths 24b, 24a, and 34, 
20 respectfully. ECM system data 31 is a subset of 
system data 18. Path 36 provides engine output 
data 38, which can be generated with each 
iteration or upon completion of a number of 
iterations as described below. Path 34 along 
25 with path 40, path 28 and path 30c, as noted 

above, deliver ECM/CM data 35, CM input data 21, 
IM/CM data 27 and CM system data 21b to CM 12. 

Path 24a, 26a and 30b direct ECM/EM data 33, 
IM/ECM/BM data 25 and EM system data 39, 
30 respectively, to EM 16 from ECM 14, IM 10 and 
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system data 18, respectively, wherein EM 16 
generates emissions data 41 through line 42. BM 
system data 39 is a subset of system data 18. 

Along with IM system data 20b described 
above, other data, ECM system data 31, CM system 
data 21b and EM system data 39, are chosen from 
system data 18. ECM system data 31 includes 
engine geometry (GEO^),RPM^. and ambient 

ten^erature (r^) . EGR valve flow area {A^,) and 
variable geometry turbine position {VGT^) may 
also be included where these features, EGR and 
variable geometry turbine (VGT) are used in the 
engine considered or controlled. Other variable 
tnay also be included depending on the engine, as 
would be understood by a person skilled in the 
art. 

Note while and are components of 
system data 18 to both IM 10 and ECM 14, they can 
be assumed or calibrated values and factored into 
10 the coinputational algorithms in both IMIO and ECM 
14 thereby avoiding actually having to provide 
these as values. This could also be said of 
GEO^, GEO^ and RPM^. however, the flexibility 
of the system is limited by setting this 
25 parameters as constants. For the purposes of 

this application, these three parameters will be 
considered provided as a subset the system data 
even if they are 'hardwired" into the modules 
that use them. 
30 CM system data 2lb and EM system data 39 
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include GEO^^ and RPM^ in both cases. 

Depending on such things as aftertreattaent in the 

case of EM 16, and the adaptability of EM 16, 

other variables may be etiployed. 

As would be landerstood by a person skilled 

in the art, GEO^g is a reference defined by a 

number of variables including such things as the 
engine cylinder diameter, number of cylinders, 

piston stroke, connecting rod length, clearance 
height, compression ratio, strokes per cycle (two 
stroke cycle or four stroke cycle) , number of 
intake and exhaust valves per cylinder, diameters 
of valves, valve lift and timing diagrams, number 
of cott5)ressors and turbines (if any - although, 
its assumed for the purposes of this disclosure 
that these are feature of the engine) , size and 
properties of cottqpressor, size and properties of 
any turbine. 

As would be understood by a person skilled 
in the art, GEO^ is a reference defined by 

variables including sizes of injector needles, 
sizes of plungers, sizes of flow channels and 
reservoirs, sizes of injector holes and number of 
injector holes. 

Output data, engine output data 38 and 
emissions data 41, includes power ['Ptntks'^ 
torque ( ) , fuel economy ( BSFC ) , exhaust 
emissions (ATO, , PM), intake air flow {Qi^) $ 
and temperature and pressure of the turbine inlet 
( ^ ' ) compressor outlet ( r„ , ) ' - 
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turbine speed (v^,) • Where used, EGR rate (Q^) 
can also be provided with engine output data 38. 
AS would be understood by a person skilled in the 
art, variables such as exhaust flow (Q^) , where 
5 important can be found. Heat release rate 

{HRR) , an intermediate parameter out of ECM 14 
originally delivered to ECM with CM/ECM data 29 
from CM 12, can also be used by the engine cycle 
model to generate engine performance parameters, 
10 such as P»,^an:d r,^, BSFC . On. , P^' T„ . and 

IM/CM data 27, IM/ECM/EM data 25 and CM/ECM 
data 29, includes rate of injection {ROD, 
determined start of injection {SOI), in the case 
15 of IM/CM data 27. In the case of IM/BCM/EM data 
25, SOI and rate of fuelling (j2jw) included. 
In the case of CM/ECM data 29, HRR and start of 
combustion {SOO is provided. These data sets, 
27, 25 and 29, form the basis of two iteration 
20 loops used to converge on a data 23 set and a 
ECM/CM data 35 set that more accurately reflects 
and refines initial input data, P^, data 19 and CM 
input data 12. ECM/CM data 35 includes the 
current iteration of , referred to as G£ , and 
25 may or may not include the current iteration of 
Q^, referred to as fi^ , depending on whether 
this is a feature of the engine used. 

ECM/EM data 33, includes intake valve 
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closing time cylinder pressure (P^^) , intake valve 
closing time cylinder temperature {T^^) as well as 
, P^, and, where a feature of the system, . 
With the flow rates of fuel, air, and EGR, 
5 pressure and teTt5)erature at the intake valve 

closure, IBM 8 module, EM 16 calculates emissions 
data 41, such as NOx emission. 

Referring to figs. 2 and- 3, small loop 98 
and big loop 118 are iteration algorithms which 
10 are used to establish initial input data of 

sufficient accuracy to allow IBM 8 to begin self 
supported operation and generate engine output 
data 38 and emission data 41, as will be 
discussed below. 
15 Referring to fig. 2, small loop 98 iteration, 
initial data, CM input data 21, CM system data 
21b and IM/CM data 27, are directed through flow 
paths 40, 30c. 28, respectfully, to CM 12 . CM 12 
calculates CM/ECM data 29 through path 32 to ECM 
20 14 where, using ECM system data 31, ECM/ CM. data 
35 is provided and directed back to CM 12 in a 
closed start-up algorithm. ECM/CM data 35 for 
each iteration, (and Q;^, where EGR is used) , 

is conipared with the previous value provided for 
25 ECM/CM data 35, 2^ (ai^ol g^' where EGR is used), 
as seen in small loop convergence step 200. 

In fig. 3, big loop 118 iteration shows the 
interaction between big loop 118 and small loop 
98. Big loop 118 includes IM system data 20b and 
30 P^i data 19 and, eventually, the current iteration 
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value of from data 23, namely, P^,. all of 

which are directed to IM 10 through paths 20, 22, 
24b. Resulting IM/CM data 27 is delivered 
through path 28 to CM 12. CM 12 uses this data, 
CM system data 21b, CM input data 21 and, 
eventually, ECM/CM data 35, delivered through 
paths 30C, 40. 34 to generate CM/ECM data 29. 
CM/ECM data 29 is delivered through path 32 to 
ECM 14. ECM 14 performs iterations of small loop 
98 to generate QH (and, where appropriate Q%,) 
until this value is converged as discussed below 
after which, ECM 14 generates P^ data 23 
delivered through path 24b to IM 10, after 
completing big loop convergence step 202. 

The data flow paths for big loop 118 are 
shown in figs. 4 and 5. Referring to fig. 4, 
initial flow path 160 is provided wherein big 
loop 118 starts from IM 10, goes through a gas 
flow loop, returning to IM 10. Subsequent big 
loop flow path 162 is shown in fig. 5. Here, IM 
10, CM 12 and ECM 14 are cycled through 
sequentially . 

In practice, referring to fig. 1, each of 
the modules 10, 12, 14 and 16, requires input 
25 data to generate results indicative of that 

aspect of the engine under consideration. The 
input data for individual modules is dependent on 
such things as the control accuracy sought, the 
demands of the user or operator and the time 
30 available for open loop control or engine. 



20 
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behavior prediction. In general, certain data is 
required in order for the controls and engine 
prediction to be practiced, as would be 
understood by a person skilled in the art. It is 
possible to supplement the controls and engine 
prediction algorithm with additional data 
depending on the needs of the user, the engine in 
question and the accuracy sought regarding the 
resulting data. 

An initial estimation is required to begin 
the subject method. The disclosure utilizes two 
iteration processes demonstrated in figs. 2 
through 5, each of which limits the impact of 
this initial estimation by limiting the 
estimation to data that is relatively easy to 
initially estimate and using big loop 118 and 
small loop 98 to correct for inaccuracies in the 
initial estimation by performing a converging 
step on this estimated values. 

IBM 8 uses system data 18, namely, geometry 
data, GEO^, GEO^ , operational data, . SOI^ . 
PW^, , RPM^ . T^, Pa^ along with initial data 
P' data 19 and CM input data 21 to generate 
output data, engine output data 38 and, if of 
5 interest, emissions data 41. CM input data 21 
may be limited to alone where no EGR is used. 
Also, where used in an engine control, the 
subject method may not, per ee, use GEOf^ , GEO^ 
but calculate output data in light of GEO^ , 
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GEO^ . That is, these geometry ^variables" may 
not be variables in an engine controller as such 
an engine would maintain the same GEO^ and GEO^ 
throughout. However, flexibility is also 
available to allow for variable valve actuation 
resulting in a variable GEO^ and variable 
compression ratio resulting in GEO^ . Each would 
be a difficult variable to map allowing IBM 8 to 
provide real time mapping for engine controls. 

T^, . and RPM^, can be measured and sent 

to a controller for control of an engine. These 
variables could be predetermined for engine 
behavio\ir prediction. Operational data, P„u • 
SOIf^ , , are operator commanded variables 

indicative of the demands required by the 
operator for the engine as translated by a 
controller or predetermined in view of a sought 
engine behaviour prediction. 

AS would be understood by a person skilled 
in the art, additional variables can be added 
depending on the engine controlled or the engine 
to be predicted. Examples include VGT^ and 
as shown in Pig. 1, reflective of an engine with 
optional components, VGT and BGR. Other variable 
would need to be included where used with the 
engine in question. 

IBM 8 relies on generation of IM/ECM/EM data 
25. CM/ECM data 29, ECM/EM data 33 by models IM 
10, CM 12 and BCM 14 to provide engine output 
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data 38 and endssions data 41. Each one of these 
data sets, 25, 29 and 33 , is dependent, however, 
on data 23, and ECM/CM data 35, generated by 
iterations perfoanned with big loop 118 and stnall 

5 loop 98. That is, IM 10 requires converged 

values for P^j data 23 on siabsequent cycles of big 
loop 118 in order to generate a data set for P^.| 
useful to provide the output data, 3 8 and 41. 
Likewise, CM requires converged values for ECM/CM 

10 data 35 on subsequent cycles of small loop 98 in 
order to generate a data set for Q^, (and fi^^ 
where EGR is used) that is useful both for 
deriving output data, 38 and 41, but also for 
supporting convergence of data 23 through big 

15 loop 118. 

ECM 14 and EM 16 need initial data in the 
method taught that is not available through the 
start-up of the method. Before ECM 14 can 
generate any useful data, it requires CM/ECM data 

20 29. Likewise, EM 16 requires ECM/EM data 33 

prior to generating emissions data 41. As such, 
there are steps within the taught method that 
must be performed before the integrated system 
will perform as a whole. Initially, IM 10 and CM 

25 12 operate with ECM in big loop 118 and small 
loop 98 within IBM 8, as described below, to 
generate the first sets of generated data, 

data 23, CM/ECM data 29, ECM/EM data 33, ECM/CM 
data 35. 
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Referring to figs. 2 and 3, big loop 118 and 
small loop 98 initiate IBM B by stepping through 
two loops. The method first performs iterations 
through small loop 98 converging on a Qi„ (and 
5 where this is used) . Once, completed, big loop 
iteration is continued, converging on a P^, . 

initiating big loop, an estimate of is 
taade, namely i^, . This allows IM/CM data 27 to be 
generated. CM 12 is then able to take initially 
10 available estimated CM input data 21, (and 
where EGR is used) and CM system data 21b along 
with IM/CM data 27 to initiate small loop 98 . CM 
12 generates CM/ECM data 29 which in turn allows 
ECM to use this data plus ECM system data 31 to 
15 provide a first iterated value of • Qfl, , 

specifically, fii. , and for each subsequent 
iteration • Again, where EGR is used, a first 
iterated value is provided, namely, ^ , and 
for each subsecjuent iteration . 
20 Convergence step 200 is performed wherein 

the relationship shown is considered. If true, 
then the convergence step is complete. More 
generally, this step involves a comparison of the 
first iteration value of (3^ , fij, , to . In 
25 general however, the iteration steps involve the 
comparison of current » fi£ ' to the previous 
value of fij,, er'- The absolute value of the 
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difference between and fi^'.^Gto. is compared to 
an error, fig, predetermined as acceptable for the 
application. An additional iteration is 
performed when 

and Qi, is set such that = when 

At this point, small loop iteration is completed 
for the first cycle of big loop 118. 
10 For the purposes of this application is 

the converged value of fiS defined as fig when 

These values are directed back through CM 12 
for subsequent iteration when \Q^-^~^\>eQ' CM 12 
15 again calculates new CM/ECM data 29, HRR and SOC 
from IM/CM data 27, SOI and iJOJ , and GEO^ and 
/?PM„g , all of which remain the same through each 

iteration of small loop 98. 

The same process is performed to determine 
20 where EGR is used. That is, the absolute 

value of the difference between and fi^'.AQ^, 
is compared to an error, Sq^, predetermined as 
acceptable for the application. An additional 
iteration is performed when 

25 !AQj=fc-!2^i>W 2 



and is set such that fi«,=fi^ when 
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For the piirposes of this application is the 
converged value of defined as fi^ when 

Upon completion of this step where Qi^ = 
(and = , where appropriate), with the 

initial iteration cycle through small loop 98 
conqQlete, iteration through big iteration loop 
continues. Here the resulting first iteration of 
P^, in this case, but generally, P^, is sent 

through convergence step 202 wherein the 
relationship shown is considered. If true, then 
the convergence step is complete. More 
generally, this involves a comparison of the 
first iteration value of P^, P^, to P^g. In 
general however, the iteration steps involve the 
comparison of current P^^, P^, to the previous 
value of P^i t P^^ . The absolute value of the 
difference between P^i and P^^ , APg^.is conqpared to 

an error, fp, predetermined as acceptable for the 
application. An additional iteration is 
performed when 

and P^ is set such that P^i^P^ when 

At this point, iteration of big loop 118 is 
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con^leted where P^, and (and where EGR is 
used) are determined allow for completion of lEM 
8 . Note that each iteration of P^i where 

5 requires a further set of iteration steps through 
small loop 98 as described above generating an 
increasingly accurate value for - 

For the purposes of this .application P^^ is 
the converged value of P^ defined as P^, when 

10 jp^.-i'^ri^^.- 

Note also, that convergence considering |AF^{ 
and |A|2ii,| can be determined by equivalent means of 
comparing iterative values of both P^i ajid 
resulting in a determination that the value has 
15 converged. These include a direct coit^arison of 
the iterative values without deteinnining the 
absolute value or the difference in the values. 
As would be understood by a person skilled in the 
art other methods are also available for 
20 assessing convergence of a sought value. 

Rxamples of suitable values for e^t , Sp 
is less than 2% of the resulting Q^, , P^ - 
The value can be much higher or much lower 
depending on the needs of the operator. 
25 The number of iteration through either big 

loop 118 or small loop 98 is, in general, 
determined by the user and the time available for 
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each iteration as balanced against the accuracy 
sought for the result for the application 
contetrvplated. This is provided by the value 
predetermined for each of Sp and . in example 

method provided. 

An advantage of the taught method is that 

the initial estimated data is limited to and 
(plus, where EGR is used, Q%,) . As these two 
values are relatively easy to estimate for engine 
controls or engine behaviour prediction, as would 
be tanderstood by a person skilled in the art, the 
efficacy of the control or accuracy of the 
prediction is determined in large part to these 
estimations and the acc\aracy of each module. 

While generally a fimction of time 
(typically expressed in degrees of crank angle 
rotation) , may also be estimated as a single 
constant pressure, generally resulting in 
accurate results after con^letion of a few 
iterations of big loop 118 through subsequent big 
loop flow path 162. In the case of engine 
contSrol, 1%, could come from the last completed 

predicted data set. 

Considering each module, IM 10 uses, along 
with P^i data 19 in a first iteration and in 
subsequent iterations, some of system data 18, 
namely, GEO^ . P„n > SOI^. PW^ . RPM^, and P^ 
to generate data 27, namely ROI and SOI . Each 
one of IM system data 20b values is provided by 
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the user as a commanded or desired or 
predetermined value. 

Turning to ECM 14, however, as well as a 
selection of data from system data 18, namely, 
GEO^, T^. Pa^ and RPM^ ^ ECM 14 also needs 
CM/ECM data 29, HRR and SOC , and IM/ECM/EM data 
25, SOI and Qjue!' These are relatively difficult 
values to estimate initially as they presume 
values dictated by an engine's upstream 
behaviour, i.e., actual SOI and not SOI^ , by way 
of example. However, CM 12 can provide such data 
along with IM 10. CM 12 uses CM input data 21, 
(and where used) . Along with this value, CM 
12 uses IM/CM data 27 as well, namely, SOI and 
ROI r and system data 18, namely, GEO^„g and 
RPM^, to provide a first iteration of HRR and 
SOC . This commences small loop 98 wherein, HRR 
and SOC provided to initiate ECM 14 as described 
above . 

Referring again to fig. 1 together with 
figs. 4 and 5, data flow through initial big loop 
flow path 160 and subsequent big loop flow path 
162, see fig. 5, illustrate the derivation of 
ECM/EM data 33 used by EM 16 along with EM system 
data 39 IM 10 generated IM/ECM/EM data 25, Qj^,^ 
and SOI, to generate emissions data 41. As none 
of IM 10, CM 12 and ECM 14 use EM 16 generated 
data, EM 16 as a module can be removed from the 
iteration process and brought back in when and 
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(and where EGR is used) are determined. 
The result of lEM 8 is the generation of 
emissions data 41 and engine output data 38, 
namely, NO,, PM , P^^, t^^, BSFC , Q^, P^, v^. 
5 T^^ Pco' Pa where EGR is used. Other 

data can be generated depending on the engine 
considered and the ECM, CM, IM or EM modules 
used. The resulting values can be compared to the 
expected values and a third iteration process can 
10 be performed converging such things as JVOj, , PM , 
P. . , r,_. , BSFC choosing different IM system data 

^ brake ' ''brake ' 

20b and different ECM system data 31 imtil the 
result sought is achieved. In effect, this is 
the iteration process used in practice where an 
15 operator will continue to demand, by way of 

example, more power continuing to adjust values 
such as , SOI^, and PW^ where these variables 
can be adjusted. For example, a driver 
commanding ever increasing power from the engine 
20 will adjust the pedal position at a rate sought 

to meet his or her needs. As would be understood 
by a person skilled in the art, the throttle 
translates into a given SOI^^ and PW^ , in 
effect, and may in some engines allow for a 
25 change in P^au ox , where used. These changing 
variables allow the engine to drive to the 
demands sought of the driver in this example. 
Similarly, a sought emissions level can be met by 
considering the results provided by EM 16 after 
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processing a first set of P„„ , SOI^ and PW^ and 
then adjusting these to meet the emissions sought 
■ while still taking into consideration the other 
demands of the pperator. in the case of 
5 controller, once the expected output data matches 
the actual output data with EGM output data 38 
and emissions data 41. the controller will 
command the engine to operate with the . SOI^ 
and PJjLf determined. Further discussion is 

10 provided below. 

Ihe engine prediction application will 
provide output data indicative of the engine 
behavior when that data meets the goals of the 
IBM operator in some respect. 
15 The following table provides an outline of 

the types of variables that can be considered to 
provide engine output data 38 and emissions data 
41 for a particular engine design: 



17ame 


Impact on; 


Injector needle diameter 


SOI , ROI , Pulse Width 


Injector tip tiole 
size/number of holes 


ROI f Qjuet 


Injector control chamber 
plunger diameter _ 


SOI , ROI ,PW,Qp^ 


Injector needle/pltinger 
travel length 


ROI ,PW,Qjuei 


Injector needle/plunger 
springs mass/stiffness 


SOI , ROI , PW , Qj,^ 


Injector sac volume 


sinr PW 


Combustion chamber 
clearance height - 
compression ratio 


P^,BSFC 
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Turbo -cnarger turome 
size 


"w -int-aW^ manifold 

pressure ( P<« ) / / Qegr 




Turbo -charger compressor 1 


Qitt f ^turb ' ^im * Qegr 




BGR configuration j 


. BSFC , m, . PM 




EGR valve size | 






Intake/exhaust valves 
timing, lift, effective 
area 


Effective Volume (,V^) , 




j^'Ff'tAy •h'T&atment devices 1 
- pipeline layout 1 


Exhaust 

tenDerattire (T^) , 


Command SOI &. PW 1 


SOI , ROI ,PW,Qfi^, NO, , PM 


Fuel rail pressure 1 


SOI , ROI ,PW .Qf^, NO, , PM 


BGR valve opening 


Q^,NO,.PM 


VGR^ coolant temperature 


EGR 

1 temperature ( ) , NO, , PM 


EGR cooler size 


T^r.Qe^>NO,.PM 


Intake/ exhaust 
restrictions 


Flow restriction of intake 
and exhaust ( Ar ) , , air 
1 to fuel ratio 
(%) , BSFC , NO, , PM 


E 1 e va t i on / envi ronment 
temperature 


\^^*,.r^.'^>BSFC.NO,.PM 


VGT^ 


Ar,Q^,7^.BSFC,N0,.PM 


Waste gate opening 


1 , A, BSFC . NO, , PM 


Valve openings of after 
treatment devices 


Distribution of exhaust 
1 flow and teinperature 



Note that for gaseous -fuelled engines, a 
pilot injection feature may be implemented within 
IM 10 with CM 12 for managing natural gas / pilot 
5 fuel combustion typical of a type of gaseous - 

fuelled direct injection high compression ratio . 
engine. The resulting CM/ECM data 29 can be 
generated to support the remainder of the method 
as would be understood by a person skilled in the 
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art. 

Referring to Pig. s, the subject invention 
is described with one embodiment for use in an 
engine control environment. Here, operator 
command 300, from an operator or driver, provides 
an engine operating command. In this example, 
pedal position, a typical from a driver, is 
translated by a look-up table, in EOT 302, into 
torque command, . This is directed with 
engine speed, RPM^. to ECU 302 from engine 305. 
see data file 304. to Supervisory Control Unit 
(SCU) 308, collectively referred to as SCU input 
data 306. SCU 308 provides, from . operating 
boundary data 310 for engine operation to ECU 
302 . These include maximum cylinder pressure 
^ pwx J ^ cgin include demanded particulate 

emissions (PAf*^), demanded NOx emissions 
( NOf^ ) , demanded tenqperatuire at turbine inlet 
(r/™*) and demanded teti5>erature at compressor 
outlet (T^). A modified torque can also be 
included where torque is to be kept within a 

certain range. 

System data is then provided to IBM 8 to 
allow it to calculate output data. Specifically, 
5 ECU 302 provides to IBM 8 a subset of input data 
312, namely an initial injector command, in this 
case, an initial start of injection {SOJl^) . an 
initial rail pressure (P,^), an initial pulse 
width (PlCd)' and and RPM^ from SCU input 
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data 306. Optionally, an initial variable 
geometry turbine position {VGT^) and initial 
effective EGR valve flow area iA^) can be 
provided where these variables reflect features 
5 of the subject engine. IBM 8 is also fed or 
stores base data 314, xiamely, GEO^ , G^EO^s 
and P^, the latter two variables provided by 
appropriate sensors. Its also possible the 
engine would be calibrated for a given and 
10 P^, in which case, an actual measurement or 
estimate of and would not be needed - 

With system data, namely input data 312, and 
base data 314, IBM 8 can then provide output 
(predicted) data 316. For the purposes of 
15 control, output data should generally include 

engine out data including r^^, an output control 
parameter which includes P^ and can include T^, 
T^t NO,, and PM and can include fuel 

consumption, usually SSFC . These values are then 
20 compared to sought torque, , and any operating 
boundary data dictated by SCU 308 in data 310 . 
In general, the demanded control parameter, 
QQjjpM ^ vhich is a f\inction of some operating 
boundary data 310, is conqpared to the generated 
25 control parameter provided in a first iteration^ 
CONP" . Two outcomes result through step 318. 
Where the relationship sought between CONP"^ , 
provided through ECU 302 as parameter data 319, 
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and CONP" and and does not satisfy the 

conditions provided for in step 318, namely: 

CONP'<CONP'^ 9 

5 where 

CONP' =f{NO„PM,P^uTco>T,) 10 

while the determined fuel consumption, BSFC , 
satisfies a minimum value to meet the conditions 
10 above {min(SSFC)) , a further iteration is 

performed through step 320 to determine a new 
predicted data set 316. Specifically, a 
convergence algorithm for deriving a multi- 
dimensional solution, such as simplex 322, is 
15 used to generate a new set of input data 324 

which is, along with and RPM^. delivered to 
IBM 8 to allow it to generate predicted data set 
316. Each n iteration through input data 324 
provides a subsequent injector command, again, in 
20 this case, start of injection (SOJ^, 

n-l,2,-.,last) , rail pressure (i^, n=l, 2,-.,last) , 
pulse width (PW^. n=l,2,>.,last) and variable 
geometry turbine position (FGT^ ,n-l,2,...,last) 
and effective EGR valve flow area 
25 {A^, n=l , 2 , last) , where appropriate . 

Note, in general, CONP*^ is a function of 
P^ as the cylinder pressure would always be 
constrained by the physical limitations of the 
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cylinder and need to be monitored. The other 
variables noted in equation 11 can be added 
depending on the engine used and the variables of 
interest. Also, while the SCU 3 08 provides the 
parameters used to determine CONP'^ . SCU 308 can 
also provide the relationship directly to ECU 302 
or IBM 8 for use in step 318. 

Eventually, step 318 will be satisfied 
resulting in step 326 through to ECU 302, where 
the last values generated for input data 324 are 
delivered to allow ECU 302 to then direct engine 
305 to operate using the last derived injector 
command, namely, a start of injection {SOI^) , a 
rail pressure {P^) and a pulse width (PW^) and 
15 a variable geometry turbine position {VGT^) and 
an effective EGR valve flow area (^^) , where 
appropriate, where each is set to the respective 
last value determined for SOP^ , P^ . PW^ and 
VGT^ and A^. where appropriate. That is. 



20 



pw^=pw!Z " 

Piens _ jylast 14 



25 In practice, such variables as the rate of 

change of pedal position are taken into 
consideration. The iterative process discussed 
above can be performed over a series of different 
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pedal positions, which should provide a rate of 
change in the pedal positions being commanded of 
the operator. 

The comparison made between the operator 
5 commands or the needs of the operator as compared 
to the conditions determined by lEM 8 as a 
possible means of meeting those needs is 
approximated within a range available for engine 
operation. That is, while the operator generally 
10 requires a given torque for his or her purposes, 
limitation of the engine physically or limitation 
set by SCU 308 for other reasons may prohibit the 
sought torque under a set of engine cominands. As 
such, ECU 302 may use, as noted above, additional 
15 boundaary data such as a modified torque able to 
approximate the operator demands and still 
operate within the physical boundaries of the 
engine and vehicle or txirbine design. 

In practice, ECU 302 can include all the 
20 components physically independent in fig. 6. In 
other words, ECU 302, SCU 308, lEM 8 and Simplex 
322 can be physically separated as suggested, but 
would likely be separate routines within the ECU 
302 itself. 

25 A further embodiment of the control strategy 

can include diagnostic capabilities. Here 
variations between predicted operating conditions 
and measured operating conditions can identify 
failures in such things as emissions sensors, BGR 

3 0 flow monitors, knock sensors and other combustion 
quality sensors, all of which could be fed, for 
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such diagnostic analysis, to ECU 302 with 
predicted data 316 and input data 324 satisfying 
step 318. 

For controls or predictive models, as would 
be understood by a person skilled in the art, the 
variables used by lEM 8 initially such as input 
data 312 or IM system data 20b can be other data 
that effectively provides the same information. 
For example while the injector comcoands of input 
data 312 include SOJ'^ . PLt > ' 
injector commands can also be reflected by an 
initial SOI^ . ROI^ , and with the same or 
similar results. Other variables can be used to 
determine and control fuel injection and 
15 combustion characteristics. 

Exanple of IBM operation 

IBM 8 integrates each of IBM's 4 modules, IM 
10, CM 12, BCM 14 and EM 16, consisting of a 
number of data files, which can be held in 
independent directories depending on their 
properties and functions. 

An example of the subject method for engine 
25 behaviour prediction is provided. The example 

does not limit the scope of the invention but is 
meant only for the purposes of illustrating an 
embodiment of the invention. For IM 10, the data 
files used include control files, f^, and P^^ 
30 files, command files, and output files. The 

control file defines SOI^ and PW^ of injection 



20 
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signals and are. therefore, a subset of IM system 
data 20b. By way of example, the following is 



typical of this file structure: 



12 



0.0 
1.000 
1.010 
1. 100 
1.110 
1.140 
1.150 
1.700 
1,701 
3.000 
3.010 
10.000 



0.0 

0.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

0.0 

0.0 



0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
2.0 
2.0 
. 0.0 
0.0 



0.0 
0.0 
50-0 
50.0 
0.0 
0.0 
9.0 
9-0 
9.0 
9.0 
0-0 
0.0 



10 



15 



Which defines the time (first column in ms) and 
commanded opening and closing timings (second 
column) where zero means closing, one means 
opening- In the example, referring to IM system 
data 20b, P,,u ^ constant, an example of the 
value for a gaseous fuelled engine being 185e5 
Pa. 

Additionally IM 10 requires, P^i data files, 
files that hold 1%, data 19. By way of example, 
backpressure file could look like: 



1440 O.S 

0.0 41.173 745.985 

0»S 42,314 751.362 

l.O 43.493 756.802 

1.5 44.711 762.306 

2.0 4S.969 767.871 

2.5 47.269 773.498 

3.0 48.611 779.186 



717.0 
717,5 
718.0 
718.5 
719.0 
719.5 



35.070 
36.005 
36.971 
37,970 
39.002 
40.069 



715.074 
720.065 
725.121 
730.240 
735.424 
740.671 
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The first line provides total nutnber of 
data points in the file and the crank angle 
interval between data points. Prom the second 
line forward, as an initial reading, the crank 
angle (CA) is provided in the first colium, and 
P»,, measured in bar, is provided in the second 
column. In this example, P^., starts from zero 
degrees CA, which corresponds to zero timing of 
the control file. A reading of 0 CA of the P^^ 
file does not mean TDC of an engine cycle, but 
means the start point of the control file. The 
injection timing corresponding to an engine cycle 
can be varied by changing the starting point of 

po 

Note as well, that P^, is, as mentioned 
above, an estimate of the backpressure in the 
system at issue. In general, this estimate can 
come from data generated on an actual engine with 
similar operational and physical specifications. 
Also, P^i could be set equal to the cylinder 
pressure determined from a previous converged 
pressure value (as is the case for Q^, see 
below) . However, where such a file is not 
available, a constant pressure, while a carude 
estimate, can be used for P^, . After iterations 
of big loop 118, P^ will provide a more 
indicative measure of backpressure during the 
course of an engine cycle on the engine at issue. 
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An example of a constant could be: 10.0e5 IcPa. 

The coramand file for IM 10, IM system data 
20b, defines the injector settings and 
operational conditions. These are operator 
driven and provide parameters of the system 
considered or system controlled. The data 
includes the variables and constants found in IM 
system data 20b and IM input data 19 as well as 
typical specifications required for the modules 
to operate and defining the boundaries of the 
iteration, for example, how long to run the 
model, the maximum time step allowed, where to 
output results . 

Also, measured or predefined values for 
UPM^ and measured values for are used. An 
example of RPM^ is 1740 and is l.oeS Pa. 
This will allow big loop 118 to commence by 
initiating IM 10. 

Also, GEO^ file could be referred to where 
the file would include, typically, information 
such as: 



Needle diameter 


[mm] 


Needle mass 


[gram] 


Needle stiffness 


[N/m] 


Plunger diameter 


[mm] 


Plunger mass 


[gram] 


Plunder stiffness 


[N/m] 


Injector hole effective 
flow area 


[mm'^2,1 


Injector hole L/D ratio 
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Needle seat angle 


[degree] 


Needle spring mass. 


[graml 


Needle spring stiffness 


[N/m] 


Plvinder spring mass 


[graml 


Plunger spring 
stiffness 


[N/m] 



Typical units are provided in the right hand 
column. 



With SOL 



P^!^ P^l^ Pet' RPMen, Pl^S 

GEOi„j , all noted above, typical output data files 
can be generated from IM 10. The main file to be 
used by lEM 8 is the formatted file, an example 
of which is provided below: 



1740.0 

-30.0 

-29.0 

-28.0 

-27.0 

-26.0 

-25.0 

-5.0 
-4.0 
-3.0 
-2.0 
-1.0 
0.0 
1.0 
2.0 
3.0 
4.0 
5.0 

65.0 
66.0 
67.0 
68.0 
$9.0 
70,0 



-4.8 



164.34 



O.OOOOOE4-00 
O.OOOOOE-l-00 
O.OOOOOE^OO 
O.OOOOOE+OO 
O.OOOOOE-fOO 
O.OOOOOE+OO 

0.00O00E4>00 
0.26492E-06 
0.18489E-05 
0.3d032E-05 
0.50160E-OS 
0.56768E-05 
0.61400E-05 
0-6S121B-05 
0.6805SE-05 
0.68804E-05 
0.69018B-05 

0,OOOQ0B4-Q0 
0.OOO0OB4-00 
O.OOOOOE-i-OO 
O.OOOOOE-t-00 
0.O000OE4-00 
0.O0OOOB<frO0 



O.OOOOOE-i-OO 
0.0OOO0E4-0O 
O.OOOOOB-i-OO 
O.0OOO0B4.00 
O.OOOOOE-i-OO 
O.0O000B-*<00 

O.OOOOOE-hOO 
O.OOOOOB-t-00 
0. OOOOOB^-OO 
O.OOOOOE+00 
O.OOOOOE-f-OO 
O.OOOOOB-fOO 
O.OOOOOE-fOO 
O.OOOOOE-fOO 
0.OO00OE4OO 
O.OOOOOB+00 
O.OOOOO&l-OO 

O.OOOOOE-I-OO 
O.OOOOOB-i-OO 
O.OOOOOE-fOO 
.COOOOOE-i-OO 
O.OOOOOE-I-OO 
0.OO000S4O0 



O.OOOOOE-l-OO 
O.OOOOOE+OO 
O.OOOOOE-I-OO 
O.OOOOOE-l-00 
O.OOOOOE-I-OO 
O.OOOOOE+00 



O.OOOOOC-t-OO 
O.OOOOOE-I-OO 
O.OOOOOE-I-OO 
O.OOOOOE+OO 
O.OOOOO&i-OO 
O.OOOOOE-fOO 
O.OOOOOB-I-OO 
O.OOOOOE-I-OO 
O.OOOOOE-I-OO 
O.OOOOOE-fOO 
O.OOOOOE-fOO 



O.OOOOOE400 
O.OOOOOE-fOO 
O.OOOOOE-I-OO 
O.OOOOOB-i-00 
O.OOOOOE-fOO 
O.OOOOOB-fOO 



O.OOOOOE-fOO 
O.OOOOOB+00 
O.OOOOOB+OO 
O.OOOOOE-fOO 
O.OOOOOE-fOO 
O.OOOOOB-fOD 

O.OOOOOB-fOO 
O.OOOOOB-fOO 
O.OOOOOE-fOO 
O.OOOOOB+QO 
0.0O00OB4-0O 
O.OOOOOE-fOO 
O.OOOOOB-fOO 
O.OOOOOE-fOO 
O.OOOOOE400 
O.OOOOOB-fOO 
O.OOOOOB-I-OO 

O.OOOOOE-fOO 
O.OOOOOB-fOO 
O.00OO0B4OO 
O.OOOOOE-fOO 
O.OOOOOE-fOO 
O.OOOOOB-fOO 



The first number being RPM^ , followed by SOI 
(measured in degree of crank angle rotation 
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relative to top dead center) , and total fuelling 
(measured in mg/shot) . The first column is 
degrees of crank angle rotation and the second 
column is ROI (measured in [kg/cal in this case) . 
Qj^ is derived from an integration of ROI data as 
would be understood by a person skilled in the 
art. The gap in data between crank angles -25 
and -5 and 5 and 65 is given for economy. The 
data file would not include these gaps. 

The result is IM 10 generating IM/CM data 27 
and IM/ECM/BM data 25. 

Having generated IM/CM data 27, CM 12 uses 
this data plus its own data file, CM input data 
21, (and if B6R is a consideration in the 
engine considered) as well as CM system data 21b, 
RPM^g and GEO^ . An exati^jle of initial CM input 
data 21 with EGR is = 1500.0 kg/hr and = 
300.0 kg/hr for natiaral gas. 

CM 12 provides CM/ECM data 29 in a file for 
use by ECM 14. An example of such a file is: 
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-20 


0 


-19 


0 


-IB 


0 


-17 


0 


-16 


0 


-IS 


0 






-5 


6.36444E-«8 


-4 




-3 


O.00O12872 


•Z 


0.00034991 


-1 


0.000618775 


0 


0.000939128 


1 


0.001366183 


2 


0.00197213 


3 


0.002863935 




0.00416602 




0.006229188 


u 






.0 




0 




. .. Q 




0 


l9B 


0 


lioo 


0 



The CA rotation relative to top dead center 
is provided in the left hand column and HRR is 
5 provided in the right hand coliiran. As would be 
understood by a person skilled in the art, SOC is 
provided with CM/ECM data 12, derived from the 
HRR data. 

Iterations through CM 12 and BCM 14 are run 
10 with small loop 98 replacing estimated Q^, 

with ei, converging to , Q^gr as taught 

above. An example of after completion of a 
small iteration loop is 1661.9 kg/hr. An example 
of the resulting data file provided to IM 10 for 
15 a subsequent big loop 162 is: 



SyxAbol 


Value 




1740.0 rpm 




1493.4 Nm 




281.5 kW 
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BSFC 


0:2177 
kg/kW h 




3.05 bar 




441.8 K 


a. 


1661.9 
kq/h 


Qjka 


51.500 
Icg/h 


SOI 


-30.0 CA 




415.0 
kg/hr 




818.0 K 




2.855 bar 



Not all of data is used by IM 10 during each 
loop as described above. The above is also 
indicative of the data generated by BCM should 
the iteration of big loop end. 

ECM 14 with HRR and SOC provides P^, an 
example of which is: 



-160 


3.02653 


-148.5 


3.03322 


-149 


3.04004 


-148.5 


3.047 


-148 


3.0541 






56r 


Z98286 


567.5 


2.99874 


568 


3.00475 


565*6 


3.01089 


569 1 


3.01716 


869^ 


""a44a7 



Where CA is provided in the first colutnn and 
P^j, in bar, is provided in the second column. 
Engine performance data, engine output data 38 , 
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ECM/EM data 33 and IM/ECM/EM data 25, is provided 
where no f\irther iterations of big loop are 
required as determined by Sp after a comparison of 
with i^"* , an exan?ple of which is : 



Sysibol 


Value 


RPM^ 


1740.0 


'"irate 


1493 • 4 JNul 




281.5 kW 


BSFC 


0,2177 




3.05 bar 




441.8 K 




1661.9 


QjUel 


51.500 


SOI 


-30.0 CA 


Qnr 


415.0 


Tn 


818-0 K 


Ph 


2.855 bar 




2.969 bar 



After iterations through big loop 118, an 
example of P,^^ from ECM 14, when P^, is set to P^, 
as described above, is: 



10 



is : 




-360.000 


1.749 


-359.500 


1.751 


-359.000 


1.756 


-358.500 


1.762 


-358. OCC 


1-770 


-357.500 


1.778 


-357.000 


1.786 


357.000 


1.784 


357. SOO 


1.773 


356.000 


1.764 


358.500 


1.757 


359.000 


X.752 


359.500 


1.749 


360.000 


1.749 



-^1 
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With the complete set of IM/ECM/EM data 25, 
ECM/EM data 16 and EM system data 39, emissions 
data 41 is generated from EM 16. An exangple of 
5 HO, based on the above is 295.3 ppm. 

A higher number of iterations through big 
loop 118 helps to refine the accuracy prior to 
fiAll generation of the ultimate output data, 
engine output data 38 and emissions data 41, from 
10 ECM 14 and EM 16. As well, the basis of this 

output data, data from IM 10 and CM 12, IM/ECM/EM 
data 25, IM/CM data 27 and CM/ECM data 29, is 
refined as well. Typically four iterations of 
big loop for an example such as the one above 
15 would be reasonable. 

CM 12 developed for directly injected 
gaseous fuel combustion is useful for the subject 
invention. The desired pressure at which gaseous 
fuels are injected into the combustion chamber of 
20 an internal combustion engine are dependent on 
to a far greater extent than is the case for 
liquid direct injection engines. Diesel fuel is 
injected at pressures that can be an order of 
magnitude higher than is the case for natural 
25 gas. As such, the range in is not as important 
to IM 10 accuracsy. Therefore, when considering a 
gaseous-fuelled engine, lEM 8 incorporating 
derived from ECM 14 after completion of big loop 
118, allows improved accuracy in IM 10 as well as 
30 iii?)roved functionality in IBM 8. Where liquid 
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fuels are used however, there is less need for 
refining the value of as is done through the 
iteration described above. As such, there has 
been less need to deal with this aspect of liquid 
5 direct injection engines. 

For the purposes of this application, the 
sub-tnodules, IM 10, CM 12, ECM 14 and EM 16 
include the prior art examples discussed above 
and TOay vary depending on the application of the 
10 user: see background discussion. 

While particular elements, embodiments and 
applications of the present invention have been 
shown and described, it will be understood, that 
the invention is not limited thereto since 
15 modifications may be made by those skilled in the 
art without departing from the scope of the 
present disclosure, particularly in light of the 
foregoing teachings. 
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Nraaenclature Table 

AUbr«viafeion Definition 

xm integrated engine model 

XPI injector module 

cooibustion module 
scat Engine cyole module 

Emiasiona module 
initial cylinder pressure 



Qt 

0° 

GEO^ 
pit 

GEO^ 

^amb 

VGT^ 

BSFC 
NO, 

PM 

Qin 

Pa 



Pco 



initial intake flow rate 
initial E6R rate 
injecstor geometry 
Iterative cylinder presaiura 
Bail pressure 

Commanded or desired start of lajection 

commanded or desired pulse width 

speed of engine 

AtiDbient pressure 
Engine geometry 

2^ient temperature 
Effective BGR valve flow area 

Variable geometry turbine position 

Brake power 
Brake torque 

Brake speoific fuel consunption 
Emissions of oocides of nitrogen 

Particulate matter emission 
intake flow rate 
cylinder pressure 
Tuzbins inlet tenqperatnre 
Turbine inlet pressure 

Con^ressor outlet: teo^rature 
GOB^xessor outlet pressure 
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HRR 
ROI 
SOJ 

Qfial 

soc 



-1 



^Qegr 

PW 

T 

Ar 

X 

■''CO 

pntax 
NOf^ 



SOU 
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csftaitien 
Turbine speed 

BGR f lov rate 

Bxhaust £lo« rate 

Beat release rata 

Rate of iajecticm 

start of injection 
Puel supply rate 

Start of coT&bixstion 

Iterative intake flow rate 

Iterative BOR flow rate 

intake valve closing time c yli nde r pressure 
intake valve dosiaa titne cylinder temperature 
previous cycle intake flow rate 
previous cycle Htai flow rate 
Intake air flow rate error 
BGR flow rate error 
Cylinder pressure error 
Pulse width. 

Intake manifold pressure 
Effective volume of cylinder 
Ekhaust teioperatttra 
Tenperature BOR 

Plow restriction of intake and eachaust 
Air to fuel ratio 

Demanded compressor outlet temperature limit 
Demanded torblne inlet tenqperature limit 
Maximum cylinder pressure 
Demanded particxilate matter emission 
Demanded emissions of oxides of nitrogen 

Commanded brake torque 

Initial commanded start of injection 
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^ ran 

vgt;„ 



CONP 



CONP 



dmd 



son 



cad 



pu 

vgt;, 

SOI„ 



PW. 



eng 



nil 

VGTp 

~ntt 

VGT^ 



Xnlt:lal commanded pulse wldtlx of Injectloa of fuel 

Znililal Qommanded rail • pressure 

Initial oooBBanded variable seometry turbine 
posit Ion 

Initial commanded effective B6R valve flow area 

Control parazoeter 

Demanded control parameter 

Iterative oononanded etart of injection 

Iterative commanded pulse widtb of injection of 

fuel 

Iterative commanded rail pressure 

Iterative eonmanded variable geometry turbine 
position. 

Iterative commanded effective EGR valve flow area 

Start of injection control command to engine 

Pulse width of injection of fuel command to enginG 

Rail preasTire command to engine 

Variable geometry turbine position command to 
engine 

Effective BGR valve flow area command to engine 
Final iterative commanded steirt of injection 
Final iterative commanded pulse widtb of injection 
of fuel 

Final iterative eoimnaaded rail pressure 

Pinal iterative commanded variable geometry 
turbine position 

Pinal iterative commanded effective IBGR valve flow 
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What is claimsd is: 

1. A method of controlling a direct injection 
internal corobuation engine, said method 
comprising 

a. calculating output data from system data, 
said system data con^jrising 
i. an injector geometry defined by a 
fuel injector in fluid communication 
with a combustion chaniber of said 



engine, 

ii. an engine geometry defined by said 
engine, 

iii . an engine operating command based on 
a desired engine output, 

iv. an initial injector command, 
V. an engine speed, 

said engine operating command being 
indicative of a commanded torque for said 
20 engine, said output data comprising 

i. an engine out value indicative of a 
brake torque delivered by said engine, 
and 

ii. an output control parameter 

25 indicative of a cylinder pressure 

value, 

b. when said output data does not satisfy a 
predetermined relationship, determining a 
subsequent injector command and 
recalculating said output data with said 
subsequent injector command, where said 
predetermined relationship compares said 



30 
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output data and said engine operating 
cotomand and a demanded control parameter, 
said demanded control parameter being 
indicative of a maximum cylinder pressure 
value , 

c. when said output data satisfies said 
predetermined relationship, cotratianding 
said injector to provide fuel to said 
combustion chamber according to a last 
injector command. 

The method of claim 1 wherein said system 
data further comprises an ambient 
ten5)erature and an ambient pressure. 

The method of claim 2 wherein said ambient 
temperature and said ambient pressure are 
measxired prior to calculating said output 
data. 

The method of claim 1 wherein said engine 
speed is measured prior to calculating said 
output data. 

The method of claim 2 wherein said engine 
operating command is converted to ^aid 
commanded torque based on an operator 
selected pedal position. 

The method of claim 1 wherein said demanded 
control parameter is determined from said 
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engine operating command. 

The method of claim 1 wherein said output 
control parameter is indicative of 
calculated emissions and said demanded 
control parameter is indicative of a 
demanded emission. 

The method of claim 7 wherein said 
calculated emissions and said demanded 
emissions are indicative of at least one of 
a NOx concentration and a particulate matter 
concentration . 

The method of claim 1 wherein said 
predetermined relationship further comprises 
comparing said output data with at least one 
previously calculated fuel consumption 
determined from said system data using said 
initial injector command or said subsequent 
injector command, said output data further 
comprising said calculated fuel consumption. 

The method of claim 9 wherein said 
predetermined relationship comprises said 
• calculated fuel consumption being less than 
said at least one previously calculated fuel 
consumption • 

The method of any one of claims 1, 2 and 9 
wherein 
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a. said initial injector coiranand comprises a 
commanded initial start of injection, a 
commanded initial pulse width and a 
commanded initial rail pressure, 

b. said siibsequent injector command conprises 
a commanded subsequent start of injection, 
a commanded subsequent pulse width and a 
commanded siibsecjuent rail pressure, and, 

c. said last injector command conqprises a, 
commanded last start of injection, a 
commanded last pulse width and a commanded 
last rail pressure. 

The method of claim 1 wherein said engine 
further comprises a variable geometry 
turbine and said system data further 
comprises an initial variable geometry 
turbine position, further comprising 

a. when said output data does not satisfy 
said predetermined relationship, 
determining a subsequent variable geometry 
turbine position, with said subsecjuent 
variable geometry turbine position and 
said subsequent injector command, 
recalculating said output data 

b, when said output data satisfies said 
predetermined relationship, commanding 
said variable geometry turbine according 
to a last variable geometry turbine 
position. 
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The method of claim 1, wherein said engine 
further comprises an EGR system and said 
system data further comprises an initial 
effective EGR valve flow area, further 
comprising, 

a. when said output data does not satisfy 
said predetermined relationship, 
determining a subsequent effective EGR 
valve flow area, with said subsequent 
effective EGR valve flow area and said 
subsequent injector command, recalculating 
said output data, 

b. when said output data satisfies said 
predetermined relationship, commanding 
said EGR system according to a last 
effective EGR valve flow area. 

The method of claim 1 wherein said output 
data is calculated using 

a. an injector module for calculating IM/CM 
data and IM/BCM data txoxa. an estimated 
initial cylinder pressure and IM system 
data, said IM system data being a subset 
of said system data, 

b. a combustion module for calculating 04/BCM 
data from an estimated initial intake flow 
rate, said IM/CM data and CM system data, 
said CM system data being a subset of said 
system data, 

c. an engine cycle module for calculating 
said output data from said CM/ECM data. 
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said IM/BCM data and said ECM system data, 
said ECM system data being a subset of 
said system data, 
wherein, prior to calculating said output 
data, a small loop between said combustion 
model and said engine cycle module is used 
to provide a converged intake flow rate and, 
a big loop between said injection module, 
said combustion module and said engine cycle 
module is used to provide a converged 
cylinder pressure, said converged cylinder 
pressure used to provide said cylinder 
pressure value. 

The method of claim 14 wherein said IM/CM 
data comprises an IM rate of injection and 
an IM start of injection. 

The method of claim 14 wherein said IM/ECM 
data comprises an IM fuel flow and an IM 
start of injection. 

The method of claim 14 wherein said CM/ECM 
data comprises a CM heat release rate and a 
CM start of combustion. 

The method of claim 14 wherein said output 
data is further calculated using an 
emissions module for calculating emissions 
data from BCM/EM data, said IM/ECM data and 
said system data, said output control 



CA 02441686 2003-09-23 



- 61 - 

parameter being indicative of calculated 
emissions and said demanded control 
parameter being indicative of a demanded 
emission. 

19. The method of claim 14 wherein said IM 
system data comprises an ambient pressure, 
said injector geometry, said engine speed 
and said injector command. 

20. The method of claim 14 wherein said CM 
system data comprises said engine speed and 
said engine geometry. 

21. The method of claim 14 wherein said ECM 
system data comprises an ambient pressure r 
an ambient temperature, said engine geometry 
and said engine speed. 

22. The method of any one of claims 1 through 21 
wherein said engine is fuelled, at least 
partially, by a gaseous fuel. 

23. The method of claim 22 wherein said gaseous 
fuel is natural gas. 

24. The method of claim 22 wherein said gaseous 
fuel conprises at least one of methane, 
hydrogen, ethane and propane. 

25. The method of claim 22 wherein said gaseous 
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fuel is hydrogen. 

A method of predicting output data generated 
by a direct injection internal combustion 
engine, said output data comprising at least 
one of engine torque, fuel consumption, 
power and cylinder pressure, said method 
comprising 

a. selecting system data, said system data 
comprising : 

i. an injector geometry indicative of a 
selected fuel injector, 

ii. an engine geometry indicative of said 
engine, 

iii. an engine speed, 

iv. an injector command, 

b. estimating initial data, said initial data 
comprising an initial cylinder pressure 
and an initial intake flow, 

c. calculating said output data using said 
system data and a converged cylinder 
pressure and a converged intake flow 
wherein said initial data used to: 

i. provide said converged cylinder 
pressure from said initial cylinder 
pressure, 

ii, provide said converged intake flow 
from said initial intake flow. 



The method of claim 26 wherein said system 
data further comprises at least one of an 
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ambient pressure and an ambient teinperature . 

The method of claim 26 wherein said output 
data further coirprises at least one of NOx 
concentration and particulate matter 
concentr at ion . 

The method of claim 26 wherein said injector 
command coinprises a commanded rail pressure, 
a commanded start of injection and a 
commanded pxilse width. 

The method of claim 26 wherein said system 
data further comprises at least one of a 
commanded variable geometry turbine position 
and an effective EGR valve flow area, 
wherein if an effective EGR valve flow area 
is included, said system data further 
comprises an initial EGR flow. 

The method of claim 27 wherein said output 
data is calculated using 

a. an injector module for calculating IM/CM 
data and IM/BCM data from said initial 
cylinder pressure and IM system data, said 
IM system data being a subset of said 
system data, 

b. a combustion module for calculating CM/BCM 
data from said initial intake flow, said 
IM/CM data and CM system data, said CM 
system data being a subset of said system 
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data, 

c. an engine cycle wodule for calculating 
said output data from said CM/ECM data, 
said iM/BOt data and SCM system data, said 
BCM system data being a subset of said 
system data, 
wherein a small loop between said combustion 
module and said engine cycle module is used 
to converge to said converged intake flow 
and, a big loop between said injection 
module, said combustion module and said 
engine cycle module is used to converge to 
said converged cylinder pressure. 

The method of claim 31 wherein said output 
data is further calculated using an 
emissions module for calculating emissions 
data from ECM/EM data, said IM/ECM data and 
EM system data, said EM system data being a 
subset of said system data, said output data 
further coTi5>rising said emissions data. 

The method of claim 32 wherein said 
emissions data ie indicative of at least one 
of a NQx concentration and a particulate 
matter concentration. 

The method of claim 31 wherein said ECM/EM 
data comprises said converged cylinder 
pressure, said converged intake flow, an 
intake valve closing time cylinder pressure 
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and an intake valve closing time cylinder 
tens)erature, said intake valve closing time 
cylinder pressure and said intake valve 
closing time cylinder . temperature provided 
by said engine cycle module. 

35. The method of claim 31 wherein said IM/CM 
data comprises an IM rate of injection and 
an IM start of injection. 

36. The method of claim 31 wherein said IM/ECM 
data comprises an IM fuel flow and an IM 
start of injection. 

15 37. The method of claim 31 wherein said CM/ECM 
data comprises a CM heat release rate and a 
CM start of combustion. 

38. The method of claim 31 wherein said IM 
20 system data comprises said injector 

geometry, said engine speed, said ambient 
pressure and said injector command. 

39. The method of claim 31 wherein said CM 

25 system data con^rises said engine speed and 

said engine geometry. 

40. The method of claim 31 wherein said ECM 
system data con5>rises said engine geometry, 

30 said ambient pressure, said ambient 

ten^jeratxire eaid said engine speed. 
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The tttethod of any one of claims 26 through 
40 wherein said engine is at least partially 
fuelled by a gaseous fuel. 

The method of claim 41 wherein said gaseous 
fuel is natural gas. 

The method of claim 41 wherein said gaseous 
fuel comprises at least one of methane, 
hydrogen, ethane and propane. 

The method of claim 41 wherein said gaseous 
fuel is hydrogen. 
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